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Abstract 

Background Neutron-deficient krypton isotopes are of particular interest due to the coexistence 
of oblate and prolate shapes in low-lying states and the transition of ground-state from one 
dominate shape to another as a function of neutron number. Moreover, the onset of large 
E2 transition strength around ^^Kr indicates the erosion of A^ = 40 sub-shell gap. 

Purpose A detailed interpretation of these phenomena in neutron-deficient Kr isotopes requires 
the use of a method going beyond a mean-field approach that permits to determine spectra 
and transition probabilities. The aim of this work is to provide a systematic calculation of 
low-lying state in the even-even 68-86j^j, isotopes and to understand the shape coexistence 
phenomenon and the onset of large collectivity around A^ = 40 from beyond relativistic 
mean-field studies. 

Method The starting point of our method is a set of relativistic mean-field-|-BCS wave functions 
generated with a constraint on triaxial deformations (/3,7). The excitation energies and 
electric multipole transition strengths of low-lying states are calculated by solving a five- 
dimensional collective Hamiltonian (5DCH) with parameters determined by the mean-field 
wave functions. To examine the role of triaxiality, a configuration mixing of both particle 
number (PN) and angular momentum (AM) projected axially deformed states is also carried 
out within the exact generator coordinate method (GCM) based on the same energy density 
functional. 

Results The energy surfaces, the excitation energies of 0^,2^*", 23" states, as well as the EO and 
E2 transition strengths are compared with the results of similar 5DCH calculations but with 
parameters determined by the non-relativistic mean-field wave functions, as well as with 
the available data. The results show a picture of oblate-triaxial-prolate shape transition. 
Coexistence of low-lying excited 0^ states is found to be a common feature in the neutron- 
deficient Kr isotopes. The underlying mechanism responsible for the shape coexistence is 
discussed. 

Conclusions The main features of the low-spin spectra and the systematics of excitation energies 
and transition strengths in the neutron-deficient Kr isotopes are reproduced very well. The 
effects of dynamic correlations and triaxiality turn out to have important influences on the 



balance between the competing oblate and prolate states. An exact treatment of configu- 
ration mixing of PN+AM projected triaxial states is highly demanded to pin down these 
effects. 
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I. INTRODUCTION 

The low-lying states of neutron-deficient even-even krypton isotopes are of particular 
interest due to the rapid structural change with the neutron number and the presence of 
multiple shape coexistence, i.e. several 0^ states with different intrinsic shapes coexist at low 
excitation energy. The irregularities observed in the ground-state bands at low spin in '^^•'^^Kr 
was firstly suggested to be explained by shape coexistence in Ref. [1]. This interpretation 
is supported by the observation of a metastable low-lying 0^ state in '^^,74^^^^ ^^ 3J g^j^^ ^j^g 
spectroscopic quadrupole moments for ''"'^•^Kr [4]. It shows clearly that the spectroscopic 
quadrupole moment of the ground-state band and that of the one based on the excited 0^ 
state have opposite signs, indicating the coexistence of a prolate ground-state with an oblate 
low-lying excited state. 

Experimental evidence indicates a rapid structural change in neutron-deficient krypton 
isotopes. From '^^Kr to ^^'^^Kr, a sudden increase of B{E2 : 0^ — ?■ 2^) value is interpreted as 
the transition from an oblate shape to a prolate one in the ground-state [5], which can also 
explain the rapid increase of charge radius [6]. From ^^Kr up to ^^Kr, the B{E2 : 0^ — t- 2^) 
value decreases smoothly, indicating the gradually decrease of collectivity towards neutron 
shell closure. The large collectivity around ^^Kr implies the erosion of A^ = 40 sub-shell gap. 

Both shell models and mean-field based methods are adopted to interpret the complex 
and rapidly changing structure of nuclei in this mass region. 

In the shell models, np-mh excitations across the A^ = 40 sub-shell are necessary to 
describe the neutron-deficient Kr isotopes, which is out of the reach of the conventional shell 
model. The complex excited VAMPIR approach using a modified G matrix with effective 
charges [7] or the Monte Carlo shell model using an effective quadrupole-plus-pairing residual 
interaction [8] has been carried out to analyze the low-lying excited states. The neutron- 
proton effective interaction and the excitation of nucleons across the A^ = 40 sub-shell into 
gg/2 orbital were claimed to be responsible for the shape mixing and the onset of large B{E2) 
values respectively. Recently, the SD-pair shell model (SDPSM) [9] with the Z = 28, A^ = 50 
core was adopted to study the collective properties in the even-even '^^-^'^Kt isotopes. The 
electromagnetic transitions among the low-lying states have been described rather well by 
a phenomenological Hamiltonian with parameters fitted to the excitation energies and by 
introducing effective charges for protons and neutrons [10]. 



In the framework of the self-consistent mean-field approaches, structural change and 
shape coexistence are described in terms of shapes and deformed shells [11]. The presence of 
multiple local minima in the energy surface associated with sizable gaps in the single-particle 
spectrum are used to interpret the shape coexistence within the mean-field approximation. 
However, these isolated intrinsic shapes could be mixed with each other under the quantum 
fluctuation. A quantitative assessment of this mixing and a detailed description of the 
spectrum require "going beyond the mean field" by means of symmetry restoration and 
taking into account fluctuations in the deformation degree of freedom. 

Shape coexistence in the neutron-deficient '^^-rsj^^j, isotopes has been studied by the mix- 
ing of particle number and angular momentum projected (GCM+PNP+IDAMP) axially 
deformed Hartree-Fock-Bogolibov (HFB) states with Skyrme force SLy6 in Ref. [12]. An 
oblate ground-state band coexisting with a prolate excited band was obtained for all light 
Kr isotopes. The failure in predicting the energy order of the oblate and prolate states could 
be due to the following reasons: 1) One is the deficiency of the SLy6 force for the fp shell 
nuclei. A prolate ground-state for '^^■'^^Kr requires a reduction of the gap at A^ = 40 and the 
separation between the /5/2 and the ^3/2 levels [12]. Besides, as pointed out in Ref. [13], the 
balance between the prolate and the oblate state can also be altered with a different surface 
energy coefficient in the Skyrme force; 2) The other is the limitation to axial states in the 
calculations. Recently, triaxiality was shown to be crucial to reproduce the transition from 
the oblate to the prolate shape from ^^Kr to '^^•'^^Kr based on the comparison between the 
collective Hamiltonian calculations with and without triaxial states for these nuclei using 
the DIS force [14]. 

In the past decades, the relativistic mean-field (RMF) theory, which relies on the basic 
ideas of effective field theory and density functional theory and is therefore referred to as 
single- reference covariant density functional theory (SR-CDFT), has achieved great success 
in the description of ground-state properties of both spherical and deformed nuclei all over 
the nuclear chart [15-18]. In particular, the energy density functional of point-coupling type 
has recently attracted more and more attention [19]. It shows a great advantage in the 
extension for nuclear low-lying excited states by implementing projection techniques [20, 21] 
and generator coordinate method (GCM) [22-25]. This types of implementation is also 
referred to as multi-reference (MR) CDFT. Since an exact MR-CDFT calculation for triaxi- 
ally deformed nuclei requires the mixing of symmetry restored triaxial states, the numerical 



computation is presently very expensive for a systematical study of nuclear low-lying states 
in medium and heavy-mass regions. As a Gaussian overlap approximation of the exact 
MR-CDFT, the model of five- dimensional collective Hamiltonian (5DCH) with parameters 
determined from the mean-field calculation is much less numerical demanding and turns out 
to be a powerful tool for the systematical studies of nuclear low-lying states [19, 26-32]. 
Within this framework, the rapid shape evolution in the neutron-rich Kr isotopes around 
A^ = 60 [33] have been studied using the PC-PKl force [34] and the results were compared 
with the similar 5DCH calculations but using the Skyrme force SLy4 [35]. 

In this paper we extend our studies to the neutron-deficient Kr isotopes around A^ = 40. 
A systematic beyond-mean-field calculation of the low-lying states will be carried out to 
examine the validity of the CDFT for nuclei with shape coexistence structure, and explore 
the mechanism for the shape coexistence, the onset of large collectivity around A^ = 40, and 
the role of triaxiality. 

The paper is arranged as follows. In Sec. II, we present a brief introduction to the 5DCH 
that is used to calculate the low-lying states in neutron-deficient Kr isotopes. In Sec. Ill, the 
energy surfaces, low-energy excitation spectra, electric monopole and quadrupole transition 
strengths are presented and discussed in comparison with the results of the 5DCH calculation 
using the non-relativistic Gogny force DIS and available data. Taking ^^Kr as an example, 
the spectra will be also compared with the exact GCM-I-PNP-I-IDAMP calculations using 
the PC-PKl force and the SLy6 force. A summary of our findings and conclusions is given 
in Sec. IV. 

II. SKETCH OF THE THEORETICAL FRAMEWORK 

The quantized 5DCH that describes the nuclear excitations of quadrupole vibration, 
rotation, and their couplings can be written in the form [26, 27, 29] 

H = Tvib + Trot + Koll J (1) 

where Kou is the collective potential, 

Koii(/3,7) = i?tot(/3,7) - AKib(/3,7) - AV,ot{f3,i), (2) 

with EtotiP,'^) the standard nuclear total energy in the mean-field calculation. The AVyib 
and AV^ot are zero-point-energy (ZPE) of vibrational and rotational motions, respectively, 
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calculated with the cranking approximation [29]. The vibrational kinetic energy reads, 
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with Jfe denoting the /c-component of the angular momentum in the body-fixed frame of a 
nucleus. It is noted that the mass parameters -B/3/3, -B/37, -B77, as well as the moments of 
inertia X^, depend on the quadrupole deformation variables [3 and 7, 

Xfc = 45fc/32sin2(7-2A;7r/3), A; = 1,2, 3. (5) 

Two additional quantities that appear in the expression for the vibrational energy: r = 
B1B2B3, and w = Bj^isB^yy — Bl , determine the volume element in the collective space. 
The corresponding eigenvalue problem is solved by means of expansion of eigenf unctions in 
terms of a complete set of basis functions that depend on the deformation variables /3 and 
7, and the Euler angles 0, 6 and ip [36]. 

The dynamics of the 5DCH is governed by the seven functions of the intrinsic deformations 
f3 and 7: the collective potential Kou, the three mass parameters: -B/3/3, B/^^, B^^, and 
the three moments of inertia X^. These functions are determined by a set of mean- field 
wavefunctions, generated by the constrainted SR-CDFT calculations on the mass quadrupole 
moments ^20 and ^22, which are related to f3 and 7 by 
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where Rq = 1.2A^^^ fm, A being the mass number. 

In the constrained SR-CDFT calculation, parity (P), x-simplex symmetry {Pe™'^''), and 
time-reversal invariance are imposed for the single-particle states. The Dirac equation is 



solved by expanding in the basis of eigenfunctions of a three-dimensional harmonic oscillator 
in Cartesian coordinate with 12 major shells, which are found to be sufficient to obtain 
reasonably convergent results for the nuclei in this mass region. If not mentioned explicitly, 
the pairing correlations are treated by using the BCS approximation with a pairing force 
separable in momentum space [37]. The numerical details of the calculations can also be 
found in Refs. [28-33, 35]. 

III. RESULTS AND DISCUSSIONS 
A. Energy surfaces 

Figure 1 displays the energy surface in the /3-7 plane for the even-even 68-86j^j, Simi- 
lar results have been obtained in the calculations with the DIS parametrization of Gogny 
force [38]. There is a triaxial minimum with (/3,7) = (0.25,42°) in ^^Kr. However, this 
minimum is soft along 7 direction and moves to the oblate side in ™Kr. Besides, a triaxial 
shoulder around (0.55, 18°) is developing in ™Kr and it becomes a local minimum in ^^Kr, 
in which, the oblate minimum becomes fragmental and is separated into two minima with 
slightly smaller and larger deformations. An abrupt structural change is shown when the 
neutron number increases from A^ = 36 to 38. The triaxial shoulder moves to the prolate 
side around (0.50,0°), which becomes the global minimum in ^^Kr. This prolate minimum 
persists in ^^Kr, but later shifts to the triaxial state around (0.40, 18°) in ^^Kr and disappears 
in the nuclei beyond N = 42. 

It is noted that the global minimum in the energy surface is shifted to the spherical 
shape at A = 40 with about 1.1 MeV lower in energy than the prolate minimum. However, 
the dynamic correlation at the beyond-mean-field level could reduce significantly the energy 
difference and leads to a largely deformed ground-state in ^^Kr. This effect will be discussed 
in detail subsequently. As the neutron number approaches A = 50, the shape of ground-state 
becomes spherical. In short, the topography of energy surfaces displays an oblate-triaxial- 
prolate-spherical shape transition picture in the neutron-deficient Kr isotopes. The energy 
surfaces by the HFB calculations using the Gogny DIS force present a similar picture [38]. 
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FIG. 1: (Color online) Energy surfaces of even-even 68-86j^j. isotopes in the /3-7 plane from a 
constrained relativistic mean-field (RMF) plus BCS calculation using the PC-PKl force. All ener- 
gies are normalized to the absolute minimum, indicated with red bullets and corresponding {/3,j) 
values. The energy difference between neighboring contour lines is 0.5 MeV. 



B. Systematics of low- lying states 

Figure 2 displays the calculated 5DCH excitation energies for 2^, 0^, and 2^ states, as 
well as the B{E2 : 2f — )■ Of) values as functions of the neutron number in the neutron- 
deficient Kr isotopes using the PC-PKl force, in comparison with available data [39]. In 
Ref. [14], similar calculations have been done using the Gogny DIS force. The Gogny DIS 
force has been adopted for a global study of the low-lying states for nuclei with proton 
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FIG. 2: (Color online) Excitation energies for 2^ , 2^, and O2 states, as well as the B{E2 : 2f — )■ 0^) 
and p^{EO : 0^ — )■ 0^) values as a function of neutron number in neutron-deficient Kr isotopes, in 
comparison with the similar calculations with the Gogny DIS force [38] and available data taken 
from Ref. [39] (squares), Ref. [4] (triangles) and Ref. [40] (circles). 

numbers ranging from Z = 10 to Z = 110 and neutron numbers A^ < 200 [38]. Therefore, 
the corresponding results are plotted as well for comparison. 

The 5DCH calculations based on the PC-PKl force give similar results to those based 
on the DIS force, both of which reproduce the systematics rather well, in particular for 
the onset of large B{E2 : 2^ -^ 0^) value around A^ = 40. Quantitatively, there are some 
differences in the prediction for the nuclei around N = 40. The PC-PKl force presents a 
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more dramatic evolution trend from A^ = 34 to A^ = 38,40. Compared with the results of 
the PC-PKl force, a lower 0^ state and a smaller B{E2 : 2^ — t- 0^^) value were obtained 
in the calculations using the DIS force. These differences are probably due to the strong 
configuration mixing in the calculation of the DIS force. 

An indicator of configuration mixing is the EO transition strength. 
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where Rq = 1.2A^^^ fm. In general, a large EO transition strength is a signature of strong 
mixing of states associated with the coexisting shapes [41]. The panel (e) of Fig. 2 shows the 
value oi p^{EO : 0^ — )> 0^^) as a function of neutron number. A reasonable agreement between 
the two 5DCH calculated results is observed in all the neutron-deficient Kr isotopes, except 
for the nuclei around A^ = 40 again. The calculation using the DIS force overestimates the 
p^{EO : 0^ -^ Ot) value in 74,76j^^ ^y ^^ f^^^^^^ ^f ^^^^^ 3^ ^Yiile the PC-PKl calculated 
results underestimates this value in ^^Kr. 

The detailed information about the structure of low-lying states is provided by the col- 
lective wave functions. Figure 3 displays the distribution of probability density pja in the 
/3-7 plane for the first two 0"*^ states and the 2f state in ™~'^®Kr, where the pja satisfies the 
following normalization condition. 



/•OO PZTT 

/ /3rf/3 / sin37ci7pj„(/3,7) = l. (8) 

^0 Jo 

As the neutron number increases, the predominate configuration of ground-state is changing 

from a weakly deformed shape (|/3| = 0.25) at A^ = 34 to a large prolate deformed one 
(|/3| = 0.50) at A^ = 38,40. The evolution of the predominate configuration shows again 
the onset of the large collectivity in the nuclei around A^ = 40. The second 0^ state is 
predominately prolate deformed at A^ = 34 and becomes a strong mixing of oblate and 
triaxial shapes at A^ = 36. The distribution of probability density pjaWil) indicates that 
the 0^ state in ™^'^'^Kt belongs to a /3-vibration state. However, the 0^ state in '^^-'^^Kr is 
a 7-vibration state and becomes a /3-vibration state again in ^^Kr. It can be understood 
from the energy surfaces shown in Fig. 1 that the '^^''^^Kr are soft under the distortion along 
7-direction. In most cases, the distribution of probability density of the 2^ state is similar 
to that of of state, except that the 2f state has already become predominantly prolate 
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FIG. 3: (Color online) Distribution of the probability density pjaiP,7) for the first two 0"*" states 



and 2+ state in ™~'^^Ki. 



deformed at A^ = 36. This picture is slightly different from that by the DIS force [14], which 
gives a predominantly oblate shape at A^ = 36. 

C. Shape coexistence and Nilsson diagrams 

Shape mixing is an important quantum concept in the study of shape coexistence in 
atomic nuclei. In order to have near-degenerate collective states with the same quantum 
numbers but different deformations, there has to be a mechanism to prevent their mixing. 
The sizable shell gaps in Nilsson diagram and large barriers between the minima in the energy 
surface is not sufficient to prevent their mixing. As discussed in Ref. [42], the intruder states 
with different parity can provide such mechanism. Therefore, shape coexistence is found 
mostly in nuclei with neutron or proton number around shell-closure and sub-shell closure 
and with the presence of intruder states [41]. 

To explore the mechanism responsible for the shape coexistence in the neutron-deficient 
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FIG. 4: (Color online) Nilsson diagram of the eigenvalues of the single-particle Hamiltonian for 
neutrons and protons in ^^Kr as a function of the axial deformation parameter (3, obtained with 
the PC-PKl force. Fermi energies are indicated with bullets. 



Kr isotopes, we plot in Fig. 4 the single-particle energies of neutrons and protons as functions 
of quadrupole deformation parameter /3. As already shown in Ref. [12], there are several 
sizable gaps located at different deformation regions in the Nilsson diagram of neutrons from 
A^ = 34 to iV = 40. 

We note that with increasing of deformation /3, the downsloping K = 1/2, 3/2 ( or 
K = 9/2, 7/2) levels from the lgg/2 orbital are diving into the Fermi sea in the prolate (or 
oblate) side and becoming intruder levels. Two large A^ = 34 gaps are located in both prolate 
and oblate sides with a similar |/3| value around 0.25, which gives rise to the mixing of a 
weakly oblate configuration with a weakly prolate one in the ground-state of ^''Kr. Since the 
intruder z/l5'9/2 orbital becomes occupied in the region with |/3| > 0.4, the mixing between 
configurations with |/3| below and above 0.4 is weak, but not zero because of pairing effect. 
Therefore, the wave function of ground-state is concentrated in the small deformed region 
(1/3 1 = 0.25), while the second O"*" state is around the large deformed region (|/3| = 0.50). 
As a consequence, the mixing between the first two 0"*" states in '^''Kr is weak, as shown in 
Fig. 3. In ^^Kr, the intruder i'lgg/2 orbital is already occupied when the deformation reaches 
to |/3| ^ 0.2, which is smaller than the \(3\ value of the state in the global energy minimum. 
As a result, there is a strong mixing between the two configurations corresponding to the 
global minimum and the second triaxial minimum respectively. There are two A^ = 38 shell 
gaps. One is located at the spherical shape, extending to the oblate side up to \/3\ ~ 0.3, 
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FIG. 5: (Color online) Low-spin spectra of '^^''^'^Kic in comparison with available data. The B{E2) 
transition strengths are given in e^ fm^. Experimental data are taken from Refs. [5, 43]. 

the other is located in the prolate side at f3 ^ 0.5. Since the prolate gap is much larger than 
the spherical-oblate one and the number of occupied intruder states is different in these two 
regions, the ground-state of '^^Kr is finally concentrated around /3 = 0.5 (cf. Fig. 3). The 
mixing between the first two 0''" states is weak, which is however not supported by the data 
ofp2(E0:0^^0^). 



D. Spectroscopy of lovi^-lying states in 
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The 5DCH calculated low-spin spectra of ™-'^^Kt are compared with available data in 
Fig. 5 and Fig. 6. The main features of the low-spin spectra are reproduced very well, in 
particular for the ground-state band and the low-lying 0^ state. The observed large E2 
transition strength from the 0^ state to the 2f state in '^'^''^^Kt was reproduced in the 5DCH 
calculations using the DIS force [4], but is underestimated in our calculations by one order 
of magnitude. It is shown in Figs. 5 and 6 that the calculated B{E2 : 0^ — )■ 2^) value 
in '^^''^^Kr by the PC-PKl force plus the separable pairing force is only half of the value of 
the neighboring nuclei, which corresponds to the weaker EO transition strength, cf. Fig.2. 
The 5DCH calculations using the SLy6 force (or using the PC-PKl force) together with the 
zero-range pairing force give much stronger EO transition with p'^{EO : 0^ — ?■ O^*") = 0.148 
(or 0.153) and larger B{E2 : 0^ — !■ 2^) value for ^^Kr. All of these results imply that the 
reproduction of large E2 transition strengths between the 2f and 0^ states requires a strong 
mixing between the first two 0^ states, as it is obtained in the 5DCH calculations based on 
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FIG. 6: (Color online) Same as Fig. 5, but for '^jejSj^j. -^^ comparison with available data taken 
fromRefs. [4, 39]. 

the DIS force. 

E. Effects of dynamic correlations and triaxiality in ^^Kr 

In this subsection, we will discuss the effects of dynamic correlations and triaxiality 
in detail by taking ^^Kr as an example. Figure 7 displays the collective potential energy 
surface, subtracted the ZPEs of vibrational and rotational motions from the total energy, 
cf. Eq.(2). It shows that the ZPEs deepen the prolate deformed minimum and reduce its 
energy difference with the spherical shape by about 0.9 MeV, which leads to a coexistence 
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FIG. 7: (Color online) (a) Collective potential energy surface of ^^Kr in the /3-7 plane [cf. Eq.(2)]. 
(b) Comparison of total energy surface and collective potential energy surface for ^^Kr. All energies 
are normalized to the spherical shape. 

picture of competing spherical and prolate minima in the energy surface. Since the prolate 
minimum is broader, the ground-state of ^^Kr is finally dominated by the large prolate 
deformed configurations in the 5DCH calculations. 

In contrast to the 5DCH results, a weak E2 transition strength B{E2 : 2f — t- Oj^) was 
obtained in the GCM+PNP+IDAMP calculation of Ref. [12] for ^^Kr using the SLy6 force 
due to the weakly oblate character of predicted ground-state band. 

Figure 8 displays the energy curves projected on particle number and different angular 
momentum, calculated using the PC-PKl force for the effective Lagrangian and a density- 
independent zero-range 6 force for the pairing channel [34]. This calculation is carried out 
using our recently developed code [23], implemented with particle number projection. The 
pfafiian method is adopted to determine the norm overlaps [44]. The numerical details 
of GCM combined with PNP+1(3)DAMP techniques based on the relativistic mean-field 
wave functions will be given in a forthcoming paper [45]. Similar results as those by the 
non-relativistic calculations in Ref. [12] are obtained. The spherical configuration is purely 
a J = configuration; hence, the projection for J = does not gain any energy. A small 
deformation of the mean-field is sufficient to introduce higher J components. The projection 
for J = gives then an energy gain, which increases rapidly to reach about 2.7(3.7) MeV 
around f3 ^ —0.15(0.20), as shown in Fig. 8. As a result, the weakly oblate deformed 
global minimum, together with a new weakly prolate deformed minimum, shows up in the 
J = projected energy surface. At larger deformation, the energy gain still increases but 
at a slower rate, which brings down the large prolate deformed minimum around /3 = 0.5. 
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FIG. 8: (Color online) (a) Projected energy curves and low-lying states from the 
GCM+PNP+IDAMP calculations for '^'^Kr. The bullets correspond to the lowest GCM solutions, 
which are placed at their average deformation, cf. [23]. (b) The energy gained from projection 
on J = 0: AE{f3) = Ej=o{/3) — E^z{P) (proj.), in comparison with the ZPEs of rotation motion 
— AVrot (rot.) and vibration motion — AVvib (vib.) evaluated with the cranking approximation. 
All the results are calculated using the relativistic mean-field wave functions by the PC-PKl plus 
zero-range pairing force. 
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FIG. 9: (Color online) Comparison of low-spin spectra from (a) the 5DCH calculations and (b) the 
GCM-|-PNP-|-1DAMP calculations using the PC-PKl plus density-independent zero-range pairing 
force for ^^Kr. 
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FIG. 10: (Color online) Same as Fig. 9, but using the SLy6 force plus density-dependent zero-range 
pairing force. The results of the GCM+PNP+IDAMP calculations are taken from Ref. [12]. 

Similar features in the rotational correction energy have been observed in the neutron- 
deficient nuclei in lead region [42, 46-48]. However, the energy gain from the projection 
for J = is not enough to make the large prolate deformed configuration lower in energy 
than the weakly oblate configuration. Therefore, a weakly deformed ground-state, instead 
of a large prolate deformed one, is also obtained in the relativistic GCM+PNP+IDAMP 
calculation. 

The main difference between our relativistic GCM+PNP+IDAMP calculation and the 
non-relativistic one is that the large deformed prolate configuration (/3 k. 0.5) here is much 
lower in energy and becomes yarst state at J = 2, in comparison with J = 4 for the non- 
relativistic one. In both calculations, the ground-state is dominated by the weakly oblate 
state, contradicting with the indication of data. As a result, the B{E2 : 2f — )■ 0^) values in 
both calculations are much smaller than the data. 

Figures 9 and 10 display the comparison of low-spin spectra from the 5DCH calcula- 
tions and GCM+PNP+IDAMP calculations using PC-PKl (with 6 pairing force) and the 
SLy6 force respectively. Similar results are found for different forces in either 5DCH or 
GCM+PNP+IDAMP calculations . In the 5DCH calculations, the energy order of the oblate 
and prolate states is correctly reproduced and the calculated electric transition strengths 
are in good agreement with the data. These results are consistent with those in the 5DCH 
calculations using the Gogny DIS force [14] and confirm the important role of triaxiality. 

We note that the dynamic correlation energies by the projections and GCM are different 
from those given by the cranking approximation. In Fig. 8, the ZPEs of rotational and 
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vibrational motions by the cranking approximation are plotted for comparison. It is shown 
that the ZPE of vibrational motion is almost constant against the deformation parameter 
P, while the ZPE of rotational motion increases smoothly with the /3. Compared with the 
energy gained from the exact AMP, the ZPE of rotational motion does not equal zero for 
the spherical state and the shell effect is not evident. As a result, the collective potential 
Vcoii, subtracted the ZPEs of rotational and vibrational motions from the total energy, does 
not present the weakly deformed oblate and prolate minima with |/3| in between 0.1 and 
0.2. From this point of view, the 5DCH calculation is easier to give a large prolate deformed 
ground-state than the exact projected GCM calculation. 

We note that in the previous studies [49, 50], the calculation of rotational correction en- 
ergy using the Gaussian overlap approximation or a so called topologically invariant Gaussian 
overlap approximation has already been performed. In particular, the rotational correction 
energy from an exact AMP calculation was obtained in Ref.[51], which demonstrated that 
the exact restoration of the rotational symmetry is fundamental for a qualitative and quan- 
titative description of the rotational energy. Therefore, an exact GCM calculation combined 
with the PNP+3DAMP techniques is required to pin down the triaxiality effect. The previ- 
ous 3DAMP studies for light nuclei [23, 25, 52] indicate that the AMP can lower the energy 
of the triaxial states. We expect that a triaxial global minimum with |/3| in between 0.2 
and 0.5 may show up in the triaxial projected energy surface with J = 0, and modify the 
conclusion drawn from the GCM calculations restricted to the axial case. 

IV. SUMMARY 

In summary, we have presented a systematical beyond RMF study of the low-lying states 
in the neutron-deficient krypton isotopes. The excitation energies and electric multipole 
transition strengths have been obtained by solving a five-dimensional collective Hamiltonian 
(5DCH) with parameters determined from the RMF calculations using the PC-PKl force. 
The results have been compared with those obtained by the similar calculations but using 
the Gogny DIS and the Skyrme SLy6 forces as well as by the exact GCM-hPNP-hlDAMP 
calculations. The results of the 5DCH calculations based on different types of forces are 
similar, except difference in the size of configuration mixing. We find a picture of oblate- 
triaxial-prolate shape transition in the neutron-deficient Kr isotopes. Coexistence of low- 
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lying excited 0^ states has been shown to be a common feature in the nuclei of this mass 
region. 

The results from the beyond mean-field calculations using the relativistic PC-PKl force, 
together with the calculated results using the SLy6 force, confirm the conclusion given in 
Ref. [14] that triaxiality is important in reproducing the energy order of the oblate and pro- 
late states. Moreover, we have analyzed the role of the intruder lgg/2 orbital in preventing a 
strong mixing between the large prolate deformed configurations with the weakly oblate de- 
formed ones, which serves as a mechanism responsible for the shape coexistence phenomenon 
in the Kr isotopes around A^ = 40. Furthermore, we have illustrated the important role of 
dynamic correlations on the onset of large collectivity in ^^Kr. 

The CDFT provides a self-consistent and universal description of nuclei all over the 
nuclear chart. This study provides an example that the CDFT can describe not only the 
nuclei with a single-configuration dominated structure, but also the nuclei with a coexistence 
structure of distinctly different shapes in low-excitation energy and the transition behavior 
from one dominate quadrupole shape to another when the effects of dynamic correlations 
and triaxiality are considered properly. Although the model of 5DCH is able to describe 
the systematics of excitation energies and transition strengths in the neutron-deficient Kr 
isotopes, the exact GCM calculations combined with the PNP+3DAMP techniques are 
highly demanded to pin down the importance of triaxiality. As the dynamic correlations 
from the ZPEs are different in detail from those by the exact projected GCM calculations, 
the conclusions drawn from the 5DCH studies might be altered. Work along this direction 
is in progress. 
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